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Abstract 



It is argued that the observed Thermodynamic Arrow of Time must arise from the 

boundary conditions of the universe. We analyse the consequences of the no boundary 

proposal, the only reasonably complete set of boundary conditions that has been put 

forward. We study perturbations of a Friedmann model containing a massive scalar field 

but our results should be independent of the details of the matter content. We find that 

gravitational wave perturbations have an amplitude that remains in the linear regime at 

all times and is roughly time symmetric about the time of maximum expansion. Thus 

gravitational wave perturbations do not give rise to an Arrow of Time. However density 
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perturbations behave very differently. They are small at one end of the universe's history, 
but grow larger and become non linear as the universe gets larger. Contrary to an earlier 
claim, the density perturbations do not get small again at the other end of the universe's 
history. They therefore give rise to a Thermodynamic Arrow of Time that points in a 
constant direction while the universe expands and contracts again. The Arrow of Time 
does not reverse at the point of maximum expansion. One has to appeal to the Weak 
Anthropic Principle to explain why we observe the Thermodynamic Arrow to agree with 
the Cosmological Arrow, the direction of time in which the universe is expanding. 
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1) Introduction. 

The laws of physics do not distinguish the future from the past direction of time. 
More precisely, the famous CPT theorem 1 says that the laws are invariant under the 
combination of charge conjugation, space inversion and time reversal. In fact effects that 
are not invariant under the combination CP are very weak, so to a good approximation, 
the laws are invariant under the time reverseal operation T alone. Despite this, there is a 
very obvious difference between the future and past directions of time in the universe we 
live in. One only has to see a film run backward to be aware of this. 

There are are several expressions of this difference. One is the so-called psychological 
arrow, our subjective sense of time, the fact that we remember events in one direction of 
time but not the other. Another is the electromagnetic arrow, the fact that the universe 
is described by retarded solutions of Maxwell's equations and not advanced ones. Both of 
these arrows can be shown to be consequences of the thermodynamic arrow, which says 
that entropy is increasing in one direction of time. It is a non trivial feature of our universe 
that it should have a well defined thermodynamic arrow which seems to point in the same 
direction everywhere we can observe. Whether the direction of the thermodynamic arrow 
is also constant in time is something we shall discuss shortly. 

There have been a number of attempts to explain why the universe should have a 
thermodynamic arrow of time at all. Why shouldn't the universe be in a state of maximum 
entropy at all times? And why should the direction of the thermodynamic arrow agree 
with that of the cosmological arrow, the direction in which the universe is expanding? 
Would the thermodynamic arrow reverse, if the universe reached a maximum radius and 
began to contract? 
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Some authors have tried to account for the arrow of time on the basis of dynamic laws. 
The discovery that CP invariance is violated in the decay of the K° meson 2 , inspired a 
number of such attempts but it is now generally recognized that CP violation can explain 
why the universe contains baryons rather than anti baryons, but it can not explain the 
arrow of time. Other authors 3 have questioned whether quantum gravity might not violate 
CPT, but no mechanism has been suggested. One would not be satisfied with an ad hoc 
CPT violation that was put in by hand. 

The lack of a dynamical explanation for the arrow of time suggests that it arises from 
boundary conditions. The view has been expressed that the boundary conditions for the 
universe are not a question for Science, but for Metaphysics or Religion. However that 
objection does not apply if there is a sense in which the universe has no boundary. We 
shall therefore investigate the origin of the arrow of time in the context of the no boundary 
proposal of Hartle & Hawking 4 . This was formulated in terms of Einsteinian gravity which 
may be only a low energy effective theory arising from some more fundamental theory such 
as superstrings. Presumably it should be possible to express a no boundary condition in 
purely string theory terms but we do not yet know how to do this. However the recent 
COBE observations 5 indicate that the perturbations that lead to the arrow of time arise 
at a time during inflation when the energy density is about 10~ 12 of the Planck density. 
In this regime, Einstein gravity should be a good approximation. 

In most currently accepted models of the early universe there is some scalar field <j> 
whose potential energy causes the universe to expand in an exponential manner for a time. 
At the end of this inflationary period, the scalar field starts to oscillate and its energy is 
supposed to heat the universe and to be transformed into thermal quanta of other fields. 
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However this thermalisation process involves an implicit assumption of the thermodynamic 
arrow of time. In order to avoid this we shall consider a universe in which the only matter 
field is a massive scalar field. This will not be a completely realistic model of the universe 
we live in because it will be effectively pressure free after the inflationary period rather 
than radiation dominated. However it has the great advantage of being a well defined 
model without hidden assumptions about the arrow of time. One would expect that the 
existence and direction of the arrow of time should not depend on the precise matter 
content of the universe. We shall therefore consider a model in which the action is given 
by the Einstein-Hilbert action 

^= Jj 1 ^ 1 ' 2 R + staL**™ 1 "" (L1) 

plus the massive scalar field action 

/* = -1/2 / d 4 x (-g) 1/2 Wd^vt + mV)- (1-2) 

JM 

In accordance with the no boundary proposal, we shall take the quantum state of the 
universe to be defined by a path integral over all compact metrics with this action. This 
means that the wave function ^[hij, </>o] for finding a three metric hij and scalar field </>o 
on a spacelike surface S is given by 

*(^-,0o) = JdW^e- 1 '^'* 1 (1.3) 

where the path integral is taken over all metrics and scalar fields on compact manifolds M 
with boundary S that induce the given values on the boundary. In general the metrics in 
the path integral will be complex rather than purely Lorentzian or purely Euclidean. 

5 



There are a number of problems in denning a path integral over all metrics, two of 
which are: 

(1) The path integral is not perturbatively renormalisable. 

(2) The Einstein Hilbert action is not bounded below. 

These difficulties may indicate that Einstein gravity is only an effective theory. Nev- 
ertheless, for the reasons given above we feel the saddle point approximation to the path 
integral should give reasonable results. We shall therefore endeavour to evaluate the path 
integral at stationary points of the action, that is at solutions of the Einstein equations. 
These solutions will be complex in general. 

The behaviour of perturbations of a Friedmann model according to the no boundary 
proposal was first investigated by Halliwell & Hawking 6 and we shall adopt their notation. 
The perturbations are expanded in hyperspherical harmonics. There are three kinds of 
harmonics. 

(f ) Two degrees of freedom in tensor harmonics. These are gauge invariant and correspond 
to gravitational waves. 

(2) Two degrees of freedom in vector harmonics. In the model in question they are pure 
gauge. 

(3) Three degrees of freedom in scalar harmonics. Two of them correspond to gauge 
degrees of freedom and one to a physical density perturbation. 

One can estimate the wave functions for the perturbation modes by considering com- 
plex metrics and scalar fields that are solutions of the Einstein equations whose only 
boundary is the surface S. When S is a small three sphere, the complex metric can be 
close to that of part of a Euclidean four sphere. In this case the wave functions for the 
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tensor and scalar modes correspond to them being in their ground state. As the three 
sphere S becomes larger, these complex metrics change continuously to become almost 
Lorentzian. They represent universes with an initial period of inflation driven by the po- 
tential energy of the scalar field. During the inflationary phase the perturbation modes 
remain in their ground states until their wave lengths become longer than the horizon size. 
The wave function of the perturbations then remains frozen until the horizon size increases 
to be more than the wave length again during the matter dominated era of expansion that 
follows the inflation. After the wave lengths of the perturbations come back within the 
horizon, they can be treated classically. 

This behaviour of the perturbations can explain the existence and direction of the 
thermodynamic arrow of time. The density perturbations when they come within the 
horizon are not in a general state but in a very special state with a small amplitude that 
is determined by the parameters of the inflationary model, in this case, the mass of the 
scalar field. The recent observations by COBE indicate this amplitude is about 10~ 5 . 
After the density perturbations come within the horizon, they will grow until they cause 
some regions to collapse as proto-galaxies and clusters. The dynamics will become highly 
non linear and chaotic and the coarse grained entropy will increase. There will be a well 
defined thermodynamic arrow of time that points in the same direction everywhere in the 
universe and agrees with the direction of time in which the universe is expanding, at least 
during this phase. 

The question then arises: If and when the universe reaches and maximum size, will the 
thermodynamic arrow reverse? Will entropy decrease and the universe become smoother 
and more homogeneous during the contracting phase? In reference [7] it was claimed that 
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the no boundary proposal implied that the thermodynamic arrow would reverse during 
the contraction. This is now recognized to be incorrect but it is instructive to consider the 
arguments that led to the mistake and see why they do not apply. The anatomy of error 
is not ruled by logic but there were three arguments which together seemed to point to 
reversal: 

(1) The no boundary proposal implied that the wave function of the universe was invariant 
under CPT. 

(2) The analogy between spacetime and the surface of the Earth suggested that if the 
North Pole were regarded as the beginning of the universe, the South Pole should be 
its end. One would expect conditions to be similar near the North and South Poles. 
Thus if the amplitude of perturbations was small at early times in the expansion, it 
should also be small at late times in the contraction. The universe would have to get 
smoother and more homogeneous as it contracted. 

(3) In studies of the Wheeler Dewitt equation on minisuperspace models 8 it was thought 
that the no boundary condition implied that \l/(a) — > 1 as the radius a — > 0. In 
the case of a Friedmann model with a massive scalar field, this seemed to imply 
that the classical solutions that corresponded to the wave function through the WKB 
approximation would bounce and be quasi-periodic. This could be true only if the 
solutions were restricted to those in which the perturbations became small again as 
the universe contracted. 

Page 9 pointed out that the first argument about the CPT invariance of the wave 
function didn't imply that the individual histories had to be CPT symmetric, just that 
if the quantum state contained a particular history, then it must also contain the CPT 
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image of that history with the same probability. Thus this argument didn't necessarily 
imply that the thermodynamic arrow reversed in the contracting phase. It would be equally 
consistent with CPT invariance for there to be histories in which the thermodynamic arrow 
to pointed forward during both the expansion and contraction, and for there to be other 
histories with equal probability in which the arrow was backward. With a relabelling of 
time and space directions and of particles and antiparticles, these two classes of histories 
would be physically identical. Both would correspond to a steady increase in entropy from 
one end of time, which can be labelled the Big Bang, to the other end, which can be 
labelled the Big Crunch. 

The second argument, about the north and south poles being similar, is really a con- 
fusion between real and imaginary time. It is true that there is no distinction between 
the positive and negative directions of time. In the Euclidean regime, the imaginary time 
direction is on the same footing as spatial dimensions. So one can reverse the direction of 
imaginary time by a rotation. Indeed, this is the basis of the proof that the no boundary 
quantum state is CPT invariant. But as noted above, this does not imply that the indi- 
vidual histories are symmetric in real time or that the Big Crunch need be similar to the 
Big Bang. 

The third argument, that the boundary condition for the Wheeler Dewitt equation 
should be VP — > 1 for small three spheres S in a homogeneous isotropic mini superspace 
model, was the one that really led to the error of suggesting that the arrow of time reversed. 
The motivation behind the adoption of this boundary condition was the idea that the 
dominant saddle point in the path integral for a very small three sphere would be a small 
part of a Euclidean four sphere. The action for this would be small. Thus the wave 
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function would be about one irrespective of the value of the value of the scalar field. 
With this boundary condition, the mini superspace Wheeler Dewitt equation gave a wave 
function that was constant or exponential for small radius, and which oscillated rapidly 
for larger radius. From the WKB approximation one could interpret the oscillations as 
corresponding to Lorentzian geometries. That fact that the oscillating region didn't extend 
to very small radius was taken to indicate that these Lorentzian geometries wouldn't 
collapse to zero radius but would bounce. Thus they would correspond to quasi-periodic 
oscillating universes. In such universes, the perturbations would have to obey a quasi- 
periodic boundary condition and be small whenever the radius of the universe was small. 
Otherwise the universe would not bounce. This would mean that the thermodynamic 
arrow would have to would reverse during the contraction phase so that the perturbations 
were small again at the next bounce. 

This boundary condition on the wave function became suspect when Laflamme 10 ' 11 
found other minisuperspace models in which a bounce was not possible. Then Page 9 
pointed out that for small three surfaces S, there was another saddle point that could make 
a significant contribution to the wave function. This was a complex metric that started 
almost like half of a Euclidean four sphere and was followed by an almost Lorentzian 
metric that expanded to a maximum radius, and then collapsed to the small three surface 
S. The long Lorentzian period would give the action of these metrics a large imaginary 
part. This would lead to a contribution to the wave function that oscillated very rapidly 
as a function of the radius of the three surface S and the value of the scalar field on it. 
Thus the boundary condition of the Wheeler Dewitt equation wouldn't be exactly \& — > 1 
as the radius tends to zero. There would also be a rapidly oscillating component of the 
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wave function. 

As before, the wave functions for perturbations about the Euclidean saddle point 
metric would be in their ground states. But there is no reason for this to be true for 
perturbations about the saddle point metric with a long Lorentzian period that expanded 
to a large radius and then contracted again. 

To find out what the wave functions for perturbations in the contracting phase are, one 
has to solve the relevant Schroedinger equation during the expansion and contraction. This 
we do in sections (3.1) and (3.2). We find that the tensor modes have wave functions that 
correspond to gravitational waves that oscillate with an adiabatically varying amplitude. 
This amplitude will depend on the radius of universe. It will be the same at the same radius 
in the expanding and contracting phases and it will be small compared to one whenever 
the wave length is less than the horizon size. Thus these gravitational wave modes will 
not become non linear and will not give rise to a thermo dynamic arrow of time. 

By contrast, scalar modes between the Compton wave length of the scalar field and the 
horizon size won't oscillate but will have power law behaviour. There are two independent 
solutions of the perturbation equations, one which grows and one which decreases with 
time. The boundary condition provided by the no boundary proposal picks out the solution 
that is a small perturbation about the Euclidean saddle point for small three spheres. It 
does not require that the perturbation about the saddle point with a long Lorentzian 
period remains small. So the no boundary proposal picks out the solution of the density 
perturbation equation that starts small but grows during the expansion and continues to 
grow during the contraction. At some point during the expansion, the amplitude will grow 
so large that the linearized treatment will break down. This however does not prevent one 
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using linear perturbation theory to draw conclusions about the thermodynamic arrow of 
time. The arrow of time is determined by when the evolution becomes non linear. The 
linear treatment and the no boundary proposal enable one to say that this will happen 
during the expansion. After that the evolution will become chaotic and the coarse grained 
entropy will increase. It will continue to increase in the contracting phase because there is 
no requirement that the perturbations become small again as the universe shrinks. Thus 
the thermodynamic arrow will not reverse. It will point the same way while the universe 
expands and contracts. 

The thermodynamic arrow will agree with the cosmo logical arrow for half the history 
of the universe, but not for the other half. So why is it that we observe them to agree? 
Why is it that entropy increases in the direction that the universe is expanding? This is 
really a situation in which one can legitimately invoke the weak anthropic principle because 
it is a question of where in the history of the universe conditions are suitable for intelligent 
life. The inflation in the early universe implies that the universe will expand for a very 
long time before it contracts again. In fact, it is so long that the stars will have all burnt 
out and the baryons will have all decayed. All that will be left in the contracting phase 
will be a mixture of electrons, positrons, neutrinos and gravitons. This is not a suitable 
basis for intelligent life. 

The conclusion of this paper is that the no boundary proposal can explain the existence 
of a well defined thermodynamic arrow of time. This arrow always points in the same 
direction. The reason we observe it to point in the same direction as the cosmological 
arrow is that conditions are suitable for intelligent life only at the low entropy end of the 
universe's history. 
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2) The Homogeneous Model. 

In this section we review the homogeneous model with metric 

ds 2 = a 2 (-N(t) 2 dt 2 + a(t) 2 dn 2 3 ) (2.1) 

where a 2 = 2/(3Trm 2 ) : N is the lapse function, a is the scale factor and <ifi| is the stan- 
dard 3-sphere metric. Expressing the scalar field as \f2^a(t> with the quadratic potential 
27T 2 a 2 m 2 (j) 2 : the Lorentzian action is 

where the dot denotes derivative with respect to Lorentzian FRW time (if not explicitely 
stated throughout the paper time derivative are Lorentzian). There are no time derivatives 
of the lapse function N in this action; it is a lagrange multiplier. Varying the action with 
respect to N leads to the constraint 

H = ^[~a 2 n 2 a + n 2 - a 4 (l - a 2 m 2 <p 2 )] = (2.3) 

where the momenta 7r a and ir^ are defined as 

a . a 3 • 

7la = ~N a = iV^ ^ ' 

and H is the Hamiltonian. This constraint is a consequence of the invariance under time 
reparametrization. Varying the action with respect to the field we obtain the reduced 
Klein-Gordon equation 

N J t (^) + 3 ^ + N 2 m 2 <t> 2 = 0, (2.5) 

This latter equation together with the Hamiltonian constraint H = 0, is sufficient to 
describe the classical dynamics. The second order equation for a can be derived from 
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these equations. In the inhomogeneous model there are also momentum constraints, but 
these are trivially satisfied in the homogeneous background. 

The quantum theory is obtained by replacing the different variables by operators. We 
will follow the Dirac method and impose the classical constraints as quantum operators. 
The Hamiltonian constraint thus becomes 



and is called the Wheeler-DeWitt equation. The solution of this equation ^o(a,(f)) is the 
wave function of the universe. There is a factor ordering ambiguity, but it is not important 
for the conclusions of our paper which rely on the classical limit. 

In this paper we investigate the predictions of the no-boundary proposal in a model 
where small inhomogeneities are taken into account. In order to impose this proposal we 
return to a path integral formulation of the wavefunction. It is very hard to calculate this 
path integral exactly. However we can have a good idea of the resulting wave function by 
using a saddlepoint approximation 



where C is a prefactor and I s ? is the Euclidean saddle-point action. In this approximation 
it is clear how to impose the proposal of Hartle and Hawking. The regularity condition is 
imposed on the (complex) saddlepoints of the path integral. The semiclassical approxima- 
tion to the path integral can then be used to estimate the wavefunction. 

One of the problems in using the semiclassical approximation in this model is that we 
cannot simply deform the complex metric into purely real Euclidean and real Lorentzian 
sections, for real arguments of the wavefunction. This could only be achieved in this 




(2.6) 



(2.7) 
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model if the time derivatives of both a and <fi vanish simultaneously on the Euclidean 
axis, 12 which is not possible as increases monotonically if the no boundary condition is 
imposed. Therefore we must solve the background equations of motion for complex values 
of time and physical variables, obtaining complex solutions which satisfy the no boundary 
proposal and have the given a and 4> on the final hypersurface. The no boundary proposal 
imposes the boundary conditions at one end of the four geometry 

a = £ = 1 ^=0 4> = 4> (2.8) 
dr dr 

thus we only have the freedom to choose the (complex) value of cj> at the origin of complex 
time r. 

Lyons 13 found that there were many contours in the complex time plane which induced 
real endpoints a and <j>. Some possibilities are obtained by choosing the initial value of <fi to 
have an imaginary part much smaller than the real part such that 4>Q m xs —(l + 2n)n/6(f)Q e 
(for integer n). In this paper we will only investigate the case n = 0. 

For small a the complex metric can effectively be considered as a small real Euclidean 
section, with </> approximately real, described by 

4> ~ 4>o an d a ~ sinm^or (2-9) 

m(po 

where r is the Euclidean time. When we consider gravitons below, it is a good approxima- 
tion to assume the following behaviour for the radius a when </>o > 1. For small a (< m</>o) 
the background is part of an Euclidean 4-sphere 

1 

a w — ■ - oo < i] E < 0. (2.10) 

mcpo cosh rj E 

The Euclidean conformal time is given by t]e = f dr/a. Although t\e has semi-infinite 
range notice that the proper distance is finite. The radius a starts at zero and increases 
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to a maximum value of l/m^o, the equator of the 4-sphere. For larger a, the saddle point 
is well approximated by de Sitter space 

a^— -± 0<r,<^--5 e (2.11) 

m(po cos 77 2 

where rj is the analytic continuation of rj E = irj. The universe is then in an inflationary 
era. In terms of comoving time: 

« O - ^ and _L e ™<M-^m 2 t 2 ( 2 .12) 

3 mcpo 

where t is the analytic continuation of r in the Lorentzian region. 
The action is given by 

^-3^( 1 -< 1 -™ 2 *' 2 > 3/2 )- < 2 ' 13 > 
For large a (>> l/m0o), the saddle point will have a large imaginary part. The wave 
function will therefore be of WKB type. After a suitable coarse graining, 14 we can associate 
the phase of the wave function to the Hamilton- Jacobi function of general relativity. When 
this is possible we will assume that the universe behaves essentially classically. The wave 
function will be associated to the family of classical Lorentzian trajectories described by 
the Hamilton-Jacobi function. 

Meanwhile the scalar field is decreasing and inflation will end at rj = tv/2 — 5 e when 
the scalar field reaches a value around unity, at which point the value of a will be a e ~ 
(l/m0o) ex P(3</>o/2)- S e is given by the implicit relation 5 e ~ exp(— 3(0o) 2 /2). For <po > 1, 
we have 5 e 1. When rj > n/2 — 5 e , the scalar field oscillates and behaves essentially as 
a pressureless fluid (i.e. dust): 

0~ 1(^)1/2 cos(mt) . (2.14) 
m a 6 
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The scale factor of the universe is then well described by 

fl~a m sinV ^ ') --6 e <7] (2.15) 

where the constants have been chosen to ensure a smooth transition between the inflation- 
ary and dust era. The universe will therefore expand to a maximum radius a m ~ mrar e ~ 
exp(9(0o) 2 /2)/m(</>o) 3 and recollapse. It will be convenient later on to redefine the origin 
of conformal time during the dust-like era by setting rjd = ij — n/2 + 35 e . The scale factor 
will then evolve as 

a « a max sin 2 — < i] d < 2n (2.16) 

Figure 1 depicts a typical classical trajectory corresponding to the no-boundary proposal. 

3) Inhomogeneous Perturbations. 

Let us now consider the behaviour of small perturbations around the the homogeneous 
model described in the previous section. We write the metric as 

9^u(t,y:) = g ltv (t) + 6g liv (t,jc). (3.1) 

The background part g^ v (t) was decribed in the previous section by the line element (2.1). 

One can decompose a general perturbation dg^ of a Robertson- Walker background 
metric into scalar (QfJ, vector ((Ptf m , (S?' e )?J tensor ((Py)^, (S°p? m , (G°f)?J 
harmonics. This classification originates from the way they transform under rotations 
of the 3-sphere. These harmonics are constructed from the scalar, vector and tensor 
eigenfunctions of the Laplacian on the 3-sphere, viz. Qf m , (S° ,e )f m and (G°j e )f m . More 
details and properties of these harmonics are given in refs. [15,16]. 
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We can expand the inhomogeneous perturbations of the metric in terms of these 
harmonics (where the index n should be thought of as a shorthand for nlm and o, e). The 
tensor perturbations are: 

Gfj are the transverse traceless tensor harmonics. The vector perturbations are: 

^ )= E^(j? r 2 J c n J S .) (3 - 3) 

where the = ST. + S 1 ^. are obtained from the transverse vector harmonics Sf. The 
scalar perturbations of the 

c W = ^ "L ( -2N*g n Q» k n Pr \ ( , 

^ Vv^V W 2a n n ij Q n + 6b n P%) ] 

where the P? = Q\J{n 2 - 1) and P% = n tj Q n /3 + Q\ l3 /{n 2 - 1) are obtained from the 
scalar harmonics Q n . We must also take into account the scalar perturbations of the scalar 
field: 

This expansion is in effect a Fourier transform adapted to the symmetry of the FRW 
background. The coefficients a ni b n ,c n ,d n , f n ,g n ,j n and /c n are functions of time, but not 
of the spatial coordinates of the three-sphere hypersur faces. Spatial information is encoded 
in the harmonics. 

In [6] the action (1.3) and (1.4) was expanded to second order around the homogeneous 
model. In appendix A, we have reproduced it with the equations of motion for the various 
Fourier coefficients. After examining the perturbed Lagrangians (A. 2) and (A. 3) we find 
that the different types of harmonics decouple from each other. Their wave functions will 
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therefore separate so we can write 



* n (a, 4>, a n ,b n ,c 



d n , fn) = Vni a i 0> a n, &n, /n)C( a > 0> C n )^ n (tt, 0, C? n ) 



(3.6) 



It is thus possible to investigate them separately. We will study the tensor and scalar 
modes in the next two subsections. For the vector modes there are only two variables c n 
and j n . The latter one however is a Lagrange multiplier and thus induces a constraint for 
the only variable left. Thus we find that the vector degrees of freedom are pure gauge and 
will only contribute to the phase of the total wave function. 

3.1) Linear Gravitons. 

Linear gravitons are the transverse and traceless part of the 3-metric and are described 
by the variables d n in the above notation. Using the background equation of motion we 
can derive the equation 17 



for the modes d n . Here the derivatives are with respect to Lorentzian conformal time 
and H = a' /a. The gravitons are decoupled from the scalar and vector-derived tensor 
harmonics and depend only on the behaviour of the background. 

We will calculate the wave function for the graviton modes using a saddle-point ap- 
proximation, assuming the background wave function (2.7) and saddle-point action (2.13). 
The tensor part of the wave function (see 3.6) can be written as 



where C = {5 2 {If t )/8d i n 5d( l ) 1 / 2 is the prefactor assuming the flat spacetime measure. 



< + 2Hd' n + (n 2 - l)d n = 0. 



(3.7) 




(3.8) 
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The Euclidean action for a mode d n calculated along an extremising path is given by 
the boundary term 



n 2 d cV v f 



lW t = (^f^+2aa'dl) * (3.9) 

1 ^E 

where r\ E is the Euclidean time, a function of the background variables a and </?o as de- 
scribed in [18]. It is possible to rewrite this action in terms of values of the field on the 
boundary d l nl d n and solutions of the classical equation p n 

d -a 2 ^-Pn- (n 2 -l)a 2 p n = (3.10) 



drj E dr) E 



evaluated on the boundary. The regularity condition for the no-boundary proposal implies 
that d n must vanish when the 3-geometry shrinks to zero and this implies that the action 
will have the form 



ir=Afi = £(^+4*)dj. (3.ii) 

2 p n a 



In regions of configuration space where the universe is Lorentzian, the appropriate analytic 
continuation of (3.11) should be taken. 

It is possible to find a good analytical approximation for p n and thus of the wave 
function using (3.8) and (3.11) and assuming that the background is described by equations 
(2.10), (2.11) and (2.15). The p n are approximately 
sinh. ji 

p n oc (cosh^ — )e nT/ £, — oo < r\ E < 0, in the Euclidean region; 

n 

sin Tj 7r 

oc (cos Ti + i )e _mT/ , < i] < 5 e , in the inflationary era; 

n 2 

cos[n(?7-37r/2 + 35 e )] sin[(r/ - 3tt/2 + 3<5 e )/2] sin[n(?? - 3tt/2 + SS e )\ \ 
X 1 cos 2 [(?7-37r/2 + 34)/2] 2ncos 3 [(?y - 3tt/2 + 3<5 e )/2] ) 

TV 

— — 5 e < rj , in the dust — like phase. 

(3.12) 



20 



Modes with nS e <C 1 are those with wavelengths much larger than the Hubble radius at 
the end of inflation. At the onset of inflation they are in their ground state and thus 
oscillate adiabatically. These modes will no longer oscillate adiabatically when they leave 
the Hubble radius during inflation. However all modes will re-enter the Hubble radius 
during the dust era when n pa tan[(r/ — n/2 + 38 e )/ 2] and start oscillating adiabatically 
again. Modes with n5 e 3> 1 oscillate adiabatically throughout the evolution. All the modes 
oscillate around the time of maximum expansion, and even if some do not have a phase 
which is exactly time symmetric, their amplitudes are. 

The variance squared of the field and its momenta for modes with nd e <C 1 around 
the time of maximum expansion are given by 

ld 2\- 1 ~ (l + 2 7 cos(2nr ? ) + 7 2 ) 

{dn) ~ 2(A* + A) ~ 2na 2 (l- 7 2 ) [6 6) 

2 A* A ^na 2 (1 - -f 2 ) 2 + 4j 2 sin 2 (2nr)) 

("dj - 2(A* + A) ~ 2 (l + 2 7 cos(2n?7)+7 2 )(l-7 2 ) ( > 

and 

i(A-A*) 47sin(2n?7) /oir , 

(d n n dn + n dn d n ) = { {A + A ^ » L _^ 2) l] (3.15) 

where 7=1 — n 2 8 2 /2. The expectation value of the Hamiltonian 

H„ = J^\7t 2 dn +4(d n n dn +n d J n )an a + d 2 J10a 2 n 2 a + 6^ (3.16) 



n ~ 2a3 



is 

Tl 

(H n ) m — at the onset of inflation 

a 

n . 

near the maximum expansion. 



(3.17) 



an 2 8 2 



This shows that modes start in their ground state before the onset of inflation and get 
excited during inflation and the dust phase. 
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A useful way to gain information about this state is to investigate the Wigner function 

Hdn^n) = ^J dAe~ 2 ™ A ip*(d n - A)ip(d n + A). (3.18) 

The Wigner function gives an idea of the phase space probability distribution of possible 
classical perturbations (once decoherence has occured). For the wave function (3.8) with 
action (3.11), it is given by 

Hdn^n) = ^r- eX P-(^^^+ ATJ^^ ~ 2i ^r^dnKn). (3.19) 

At the onset of inflation the Wigner function is a round Gaussian (factoring out the mode 
number and the radius of the universe). A mode with n < tan(7r/2 — 5 e ) will go outside 
the Hubble radius and have frozen amplitude and the Wigner function will then become 
an ellipse elongated in the momentum direction. When the mode comes back within the 
Hubble radius it starts rotating with period 2n/n in phase space. This behaviour lasts 
until n m tanry in the recontracting phase. The parameter characterizing the eccentricity 
of this ellipse is called the squeezing and has been studied by Grishchuk & Sidorov 19 . 

Typical classical perturbations d^ resulting from the above Wigner function are small 
at the onset of inflation. Their amplitudes get frozen when they leave the Hubble radius. 
During this stage their energies increase. The perturbations will start oscillating again 
with amplitude proportional to a -1 when they come back within the Hubble radius in the 
dust phase. They behave like 

d simmy + e) 
dn ~ anV*5 e (3 - 20) 

where e is an unimportant phase depending on the details of the matching of the p n 

functions in (3.12). Around the time of maximum expansion the amplitude of the graviton 
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modes is symmetric and thus their arrow of time agrees with the cosmological one. Figure 
2 depicts a typical classical evolution of a linear graviton. 

3.2) Linear Scalar Perturbations. 

(a) Quantum Mechanics of the Physical Degree of Freedom 

We have seen that gravitons are adiabatic near the time of maximum expansion so 
that their amplitude is time symmetric with respect to that point. This is not special 
to gravitons as the electromagnetic field, massless or conformally coupled scalar fields will 
also be adiabatic. In this section we will show however that perturbations of massive scalar 
field will not behave adiabatically at the time of maximum expansion. 

From the expansion (3.4) and (3.5) we see that there are five scalar degrees of freedom 
described by the time-dependent coefficients a n ,b n , f n , k n and g n . However the latter two 
appear as Lagrange multipliers in the Lagrangians (A. 2), (A. 3) and induce two constraints 
so overall there is only one true scalar degree of freedom. Without the presence of the 
scalar field the scalar degrees of freedom would also be pure gauge. Care should be taken 
in the treatment of the scalar perturbations in order to avoid gauge dependent results. Let 
us first find the real degree of freedom. 

Variations of the action with respect to the Lagrange multipliers N, g n and k n result 
in the Hamiltonian, linear Hamiltonian and momentum constraints. In Dirac quantization, 
which we follow here, these constraints are imposed as constraints on the quantum state. 
The wave function therefore depends only on a linear combination of the coffecients a n , b n 
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and f n . The momentum constraints ensures that the wavefunction is invariant under 
diffeomorphisms of the spatial three- surfaces. The Hamiltonian and linear Hamiltonian 
constraints ensure time reparametrization invariance of the wave function. 

Shirai and Wada 20 give an explicit form for the wave function which automatically 
satisifies the momentum constraints. These are solved by making the judicious change of 
variables 

n V > (321) 

= - 3^2b n f n . 

n 

where a = In a. Once this transformation has been performed, the momentum constraints 
imply that the wave function is independent of the linear combination a n — b n . In terms 
of the two degrees of freedom left, the linear Hamiltonian constraint becomes 

7i>7r /n - n a n Sn + e 6a m 2 0/ n + K n s n = (3.22) 

where s n = a n + b n and K n = i[(n 2 -4)7r 2 - (n 2 + 5)ir 2 - (n 2 -4)e 6a m 2 (/) 2 ]. The remaining 
gauge degree of freedom can be eliminated by solving the linear Hamiltonian constraint 
using the canonical transformation 

(y n \ = (K n e 6a m 2 <A (s n \ 

\Z n ) \1^4> 7T Q J \fn J 

(3.23) 

(*s n \ = ( K n nA(n yn -^\ 
\n f J [e^m 2 ^ ttJ ^ n Zn J 

where E = —K n Tv a + e 6a m 2 07r^>. The linear Hamiltonian constraint then implies that, 
imposed as a quantum constraint, 
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(3.24) 



and so \1/ is independent of y n . Therefore the true degree of freedom has been isolated 
the wave function is found to depend only on the single physical variable 



Z n = ^4>S n + K a f n = a 2 ((p's n - Hf n ) 



(3.25) 



and on the background variables a and 4> (in the rest of the paper we will drop the tilde 
on a and </>). The expression for the Hamiltonian for the modes z n is rather complicated 
and is shown only in Appendix A. 

We can find the the wave function for the scalar perturbations in terms of the real 
degree of freedom by using the semiclassical approximation to the path integral expression 
for the wave function as in the graviton case 



^ s (a, 0, z n ) ~ C(a, (p) exp(-J 



(3.26) 



The Euclidean action of the saddlepoint contribution to the path integral is a boundary 
term (since the action is quadratic) given by 



I c n l = (Mz n z' n -Nz 2 Xf 



(3.27) 



for 



M = 
N = 



(n 2 - 4) 



2[(n 2 -4)a' 2 + 3a V 
1 



K n (2a 4 -3a 6 mV+3 (n2 " 1) ~ 4 -' 2 



aV) + a 12 mV + 3a 



(n 2 - 4) 



b'a' 



4MUa 3 
U = K n aa' + a s m 2 (j)(j)' 

and the derivatives here are with respect to Euclidean conformal time. 

It is difficult to find solutions of the equation for z n . It is easier to return to the 

original variables and pick a particular gauge. In order to study the scalar perturbations 
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we shall choose the gauge b n = k n = 0, which is known as the longitudinal gauge. Once 
the result has been obtained in this gauge it will be easy to recast it in terms of the true 
degree of freedom z n and therefore in a gauge invariant way. Alternatively, we could use 
the gauge invariant variables of Bardeen. 21 Their relationship with the formalism used here 
is described in appendix B. 

In the b n = k n = gauge we have the equations of motion (in Lorentzian time) 



< + 3Ha' n + (3m 2 <j) 2 a 2 - 2)a n = 3(m 2 <j)a 2 f n - <f>' ' f n ) (3.28) 

fl + mf n + (n 2 - 1 + m 2 a 2 )f n = 2m 2 <f)a 2 a n - A(j)'a' n (3.29) 
and the constraints 

a' n + Ha n = -3(j> , f n (3.30) 

a n (n 2 -A- 3<p' 2 ) = Z<j>' f n + 3m 2 <Pa 2 f n + 9H<P' f n . (3.31) 



Equations (3.28-30) are just (A. 7), (A. 11), (A. 8), noting that g n = — a n in this gauge. The 
last equation follows from (A. 14), (3.30) and the background constraint. These equations 
are not independent, the first one can be obtained by taking a derivative of the first con- 
straint and using the second equation and the background equation of motion. Equations 
(3.28), (3.30) and (3.31) can be combined to give the decoupled equation of motion for a n : 

< + 2(H - ^tK + (2H' - 2H^- + n 2 + 3)a n = (3.32) 
(p (p 

This equation is useful in the inflationary era where 0' ^ 0. It is also useful in the limit 
where the curvature of the 3-space can be neglected as we can solve it explicitly in either 
the adiabatic or non-adiabatic regime (see [22]). Once we have a solution for a n , we can 
also find f n using the constraint equations (3.30) or (3.31), and therefore the real degree of 
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freedom z n . In the region near the maximum expansion it is much harder to solve (3.32) 
and we return to (3.28-3.31). 

(b) No Boundary Proposal Mode Function 

Let us now construct the solutions of (3.28-3.32) selected by the no boundary proposal. 
We focus only on modes which go outside the Hubble radius during inflation. These are the 
ones which get excited by the varying gravitational field. The very high frequency modes 
remain adiabatic throughout the history of the universe, so their arrows of time will agree 
with the cosmological one. As in the graviton case we divide the background saddle-point 
4-geometry into an approximately Euclidean section, followed by an inflationary one which 
finally turns into dust. We have however to take into account the detailed behaviour of 
the background scalar field cj> as it couples directly to the perturbations. We first find the 
regular Euclidean solutions and match them up to the ones in the inflationary phase. This 
can be done by analytic continuation. In the inflationary era the modes oscillate for a 
while until they leave the Hubble radius. At that point we match them to nonadiabatic 
solutions. Finally, the inflationary era comes to an end when <fi becomes small and starts 
oscillating, behaving like a dust background. At this point we match on the solutions 
for the dustlike phase. It turns out that for the Euclidean and inflationary solutions the 
right hand terms in (3.29) are negligible. We can solve for the scalar field modes f n and 
calculate a n from an integral version of the constraint (3.30) and check that this agrees 
with the approximate solutions of (3.32). If these terms were negligible during the whole 
of the dust era the modes would oscilllate adiabatically around the maximum expansion as 
in the graviton case. However we show that these terms do contribute to a monotonically 
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increasing amplitude of the scalar field perturbations around maximum expansion. 

The no boundary proposal requires that the matter fields in the path integral be 
regular, so in the semiclassical approximation we look for solutions to the Euclidean per- 
turbation equations which are regular as r — > 0. The regularity condition requires that 
/„ and a n vanish as r — > 0. For n ^> 1, the dominant terms of equation (3.32) are the 
second derivative of a n and — n 2 times a n and one can construct a WKB solution. The 
approximate Euclidean solution selected by the no boundary proposal is 

a n **A^e nr >*, f n ^-^— (3.33) 
a 3 a 

for some complex constant A. Here, the conformal time tje = corresponds to the juncture 

of Euclidean and Lorentzian spacetimes. Continuing the regular Euclidean solution into 

the Lorentzian section, taking r\ E — > in, gives 

1 • An e inr] 

a n » -imAe mr > f n » (3.34) 

3 3 a 

where we have used <j>'/a = im/3 during inflation (dash now denotes Lorentzian time 

derivative). The analytical continuation holds into the inflationary era as long as the 

wavelength is smaller than the Hubble radius, i.e. n 3> H. By this time inflation has 

begun and we can match onto the inflationary solutions. When the modes move outside 

the Hubble radius the modes a n and f n stop oscillating. They both have decaying and 

growing modes (the latter would be constant in the limit of exact de Sitter space). As the 

universe inflates only the slowly growing mode remains 22 so that 

On ~ f n ~ (3.35) 

where D = \m%Ae mr]H (4>% + ~~ ) is a constant depending on the detailed matching of 
the modes when they cross the Hubble radius at the time n H . This solution is valid until 
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the background scalar field decrease to <f> ~ 1. Figure 3 depicts the behaviour of a n during 
inflation and the beginning of the dust phase. 

Eventually inflation ends and the background scalar field begins to oscillate. We 
expect that the the background will behave effectively as a dust-filled universe (see equation 
(2.16)) for perturbation modes with physical wavelengths much larger than the scalar field 
compton wavelength (n <C ma) since the pressure of the oscillating scalar field averages to 
zero over that wavelength scale. Therefore the metric perturbations will behave like those 
of a pure dust universe (see, e.g. [22]). This is indeed what is found below. 

During inflation the Hubble radius H~ l is roughly constant but as the universe evolves 
in the dust era the Hubble radius starts growing. When it becomes larger than the compton 
wavelength 1/ma, the dominant term in (3.29) is m 2 a 2 . The perturbation of the scalar 
field will start oscillating again. In this early stage of the dust era when the curvature of 
the 3-surface is negligible it can be shown that the f n oscillate exactly in phase with 0' as 
follows: 



This will remain true in later stages of the dust era as long as n < ma e . This condition 
ensures that the phase of f n obtained by integrating (3.29) does not differ appreciably from 
that of 4>'. Using (3.36) together with (3.30) we can establish that the metric perturbation 
a n , time averaged over one oscillation period of n/m, is growing. The small oscillations 
around this average arise because the background energy momentum tensor is not exactly 
that of dust but that of an oscillating scalar field. The averaged gravitational perturbation 
can be calculated by taking the derivative of the averaged version of (3.30) to obtain 




(3.36) 



29 



the differential equation 

af + SHa£' - 2a£ = 0. (3.37) 



The general solution is a linear combination of the solutions 



anti _ S111 ^ (o oo\ 

°» ~(l-COS % )3 (3 - 38) 



and 



sym _ 2 sin 2 r] d - 6(rj d - tt) sin rj d - 8 cos rj d + 8 
a ™ (l-cosr^ • 

The conformal time is defined with the new origin at the beginning of the dust phase 

(i]d ~ 0) . These solutions are antisymmetric and symmetric with respect to the maximum 

of expansion (t] d = tt) and are the same solutions found for perturbations in a pressureless 

perfect fluid universe, as expected. Both solutions diverge like rj^ 5 in the beginning of the 

dust era as i]d — > 0. There is however a regular solution, given by a^ eg := a^ ymm — 67ra^ ntl , 

which approaches a constant in this limit. At the end of inflation the a n picked out by 

the no boundary proposal are small as seen from (3.35). Therefore the regular solution is 

the one selected by the no boundary proposal and this is asymmetric in the dust era: the 

perturbation amplitude steadily increases with time. Matching the solutions for the dust 

era to (3.35) shows that during the dust era 

a n w ZK eg (3-40) 

We can now use (3.36) to see that f n is oscillating with monotonically increasing amplitude 
throughout the dust era: 

fn ~ ~Tx — — — vi^d - 6sin?7 d + 2r] d cosr] d } (3.41) 
(l-cos?? d ) 2 
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With these solutions we can construct the wave function (3.26). When the background 
saddle-point is approximatively Lorentzian, the no-boundary wave function for the scalar 
perturbation is 



where M is given in (3.27) and \x n is the modefunction for z n . It is a solution of the 
equation of motion for z n picked out by the no-boundary proposal. It is explicitly given 
by the function a n and f n using (3.25) with z n replaced by /i n . From the solution of a n 
and f n we can see that it is clearly asymmetric about the time of maximum expansion. 
Considering points placed symmetrically about the maximum of expansion, the background 
will be the same at both points so that the asymmetry in the modefunction manifests itself 
as an asymmetry in the wavefunction. The variance of z n is proportional to the modulus 
of \x n and is therefore asymmetric with respect to the time of maximum expansion. We 
therefore conclude that the wavefunction predicts the continuing growth of low frequency 
scalar perturbations even when the universe begins to recollapse. 

This result alone provides a time asymmetry so long as the modes stay in a regime 
where they can be treated in a linear approximation. However, most modes will also enter 
a nonlinear regime well before the maximum expansion occurs. When this occurs the 
interaction terms in the Hamiltonian will become important and hence the coarse-grained 
entropy will increase throughout the evolution. 

Considering the stress tensor in the gauge-invariant formalism (see e.g. [22 ]), we can 
show that the density contrast is 




(3.42) 



bp 



2a 



[(n 2 -4)a 



9W/n] 



(3.43) 



3a 



n 



P 



rn 
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Modes will cross the horizon (H ~ n) when rjd ~ 1/n, and the recent COBE results 5 
indicate that the density contrast at this time is of order 10~ 5 . Using equations (3.40) and 
(3.41) in (3.43), we find that the constant D is of order 10~ 5 . At later times, only the first 
term in (3.43) is important and we find that the density contrast behaves like 

6 -j « 10" 7 n 2 ^ (3.44) 

Consequently, when the density contrast is of order unity we expect nonlinearity to be the 
dominant feature and this occurs for rf d > 10 7 /n 2 . Modes with n > 1000 will therefore 
enter a nonlinear phase before they reach the maximum and the coarse-grained entropy 
for these modes will grow. 

4) Conclusion. 

In this paper we have investigated the consequences of the no-boundary proposal 
for the arrow of time. In particular we have investigated the behaviour of small metric 
and matter perturbations around a homogeneous isotropic background. The no-boundary 
proposal predicts classical evolution with an inflationary era followed by a dustlike era. We 
found that perturbation modes are in their ground state at the beginning of the inflationary 
era. This can be interpreted as a statement that the universe is born in a low entropy state. 
Modes which leave the Hubble radius during inflation become excited then subsequently 
evolve in various ways in the dustlike era. 

We find that gravitons oscillate adiabatically for most of the dustlike era and conse- 
quently the amplitude of their oscillations is time symmetric with respect to the point of 
maximum scale factor. However, looking at the physical scalar degrees of freedom we find 
that those which have been excited by superadiabatic amplification during inflation have 
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a time asymmetric evolution with respect to the maximum. In particular, the variance of 
the scalar modes predicted by the wavefunction is different at the same value of the scale 
factor before and after the maximum. 

Thus we find that the wavefunction of the universe distinguishes between symmet- 
rically placed points on either side of maximum volume. The expanding phase has a 
smaller amplitude of the variance in the low frequency scalar modes than does the cor- 
responding point during the collapsing phase. In other words, the thermodynamic arrow 
coincides with the cosmological arrow before the maximum, but points in the opposite 
direction after the maximum. This is true for all the lowest frequency modes, so that they 
induce a well-defined thermodynamic arrow of time. Amongst the modes which display 
this nonadiabatic behaviour, higher frequency modes will enter a nonlinear regime dur- 
ing the expansion and consequently produce a growing coarse-grained entropy throughout 
expansion and recontraction, and hence also create a thermodynamic arrow of time. 
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Appendix A) Action and field equations. 

In this appendix we reproduce the action and field equations of the perturbed FRW 
model driven by a massive minimally coupled scalar field from ref.[6]. The homogeneous 
part of the Einstein-Hilbert Lagrangian is 

io-^o^-i-^^V] (Al) 
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The second order perturbation of the Einstein Hilbert Lagrangian is 



= V2aN { ^- u -)al+ [n \ 7) ^ _ 4 i ) ) 6^-2(n 2 -4)4-(n 2 + l)^ + ^(n 2 -4)a n 6 r 



1, o 5, 



+ ^n[(n 2 - 4)6 n + (n 2 + l/2)a n ] + 



3(n 2 - 1) 



k 2 n + (n 2 - 4) j; 



iVo 



(™ 2 " 4) , 



a . a 



+1/2 -^r- { - d 2 + f-= -{b 2 n + (n 2 - 4)c 2 + d 2 n + g n [2-d n + ^(3a n - g n )\ 



(n 2 - 1) 



a ar 



+ - 



( n 2 _ 4 \ 

- 2a n d n + 8—^ -rb n b n + 8(n - A)c n c n + 8d n d r 

(n z — 1) 



+ - 



--ai+6 



(n 2 - 4) 
(n 2 - 1) 



b 2 n + 6(n 2 -4)c 2 +6d 2 ] 



1 

+ - 

a 



^7 f • (™ 2 - 4 ); a 1 rt/ 2 ^- • 

-k n \ ~a n - — ~ — b n + -g n \ - 2(n - 4)c n j n 

3 I [n z — 1 a J 



(n 2 - 1) 



(A2) 



The perturbation of the matter Lagrangian gives: 



L n m = \/2N aH -j^(f 2 + 6a n f n <P) - m 2 (f 2 + 6a n / n 0) - ^(n 2 - l)f 2 + ^g 2 



2/ j-2 



1 ' 2 lW 2 . <F 2 



+ - 



— m 



2 ±2 



fl! _ 4 (n 2 -4) 



(n 2 - 1) 



b 2 n - 4(n 2 - 4)c 2 - 4d 2 ] 



~9r, 



2m 2 f n (j) + 3m 2 a n (j) 2 + 2%f + ^ 



N 2 N?: J aN, 



2 Ar 2 knfnS 



(A3) 



The field equations necessary to calculate the saddle point approximation are given 



below. From (A.l) we find the equations obeyed by the homogeneous background fields. 



The homogeneous scalar field </?o obeys 



N °it 



1 d(po~\ da dcpo 



-N dt 



+ 3 



adt dt 



+ Nom 2 (po = quadratic terms, 



(A4) 
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and the scale factor a obeys 



N °7t 



1 da 

N Q a dt 



N 2 3 d 2 N 2 
+ 3<^o t - t, ( o + r + Nlm 2 (f) 2 ) = quadratic terms, (A5) 



a 2 2 V a< 



or 



The background variables a, <^o and their momenta are subject to the constraint 
0o 1 



a 



+ o + m = quadratic terms. 



a 2 iV 2 iV, 2 a 2 



(A6) 



Let us now turn to the equation of motion of the small inhomogeneities. Variations 
with respect to a n , b n , c n , d n and f n give the following second order field equations: 



,r d 
N °dt 



a 3 da r 



1 



v /f + ^(n 2 - 4)N 2 a(a n + b n ) + 3a 3 (0 o /n - N 2 m 2 ip f n ) 



N 2 [3a 3 m 2 ip 2 - ^(n 2 + 2)a]g n + a 2 dg n - ^iV -| 



ci 2 k n 



a 3 db n 
No~~d7 



^(n 2 - l)N 2 a(a n + b n ) = ^N 2 (n 2 - l)ag n + -N ^ 



N Q 

d [a 2 k r 



N 



o 



(A.7) 
(A8) 



No 



Nr 



d_ 

~dt 

d_ 

] Jt 



' a 3 dc n ' 


d 




IN dt - 


~ dt 


[ No \ 



a 3 dd r 
PJo~~dt J 



+ (n 2 - l)N 2 ad n = 0, 



(A9) 



(A10) 



and 



Nn 



dt 



a 3 df n 
N dt 



+ 3a 3 <p d n + Nq 2 [m 2 a 3 + (n 2 - l)a]f n 



3/ o;\r2 2 , • • VO^nx 

= a {—2N Q m (p g n + Vo9n ) 



(All) 



Variations with respect to k n ,j n and g n lead to the constraints 



■ ^ ~ 4 ) L , • ^ 

(n z — 1) a 



a(n 2 -l)' 
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(A12) 



(A13) 



and 



a z ■ an 9M 2 1 

Sa n (<p 2 - ^) + 2(^ /n - ^) + iV 2 m 2 (2/ n ^ + 3a n ^) - -^[(n 2 - 4)6 n + (n 2 + -)a n ] 



"3^2" 



a 



0/ 



(A14) 



We also give the perturbation Hamiltonian in terms of the real degrees of freedom: 



H%(z n , 7T Zn ) = Ayr 2 + Bz n TV Zn + Cz, 



with 



. 1, 2 3a 3 <p 2 , 
A = -(aa 2 + y ) 
2 V (n 2 -A)' 



B = — 



C 



1 

1 



2U 2 



in 2 — 11 

K(2a - 3a 3 m 2 </> 2 - 3 7-= TraV) + a 9 m 4 (/) 2 - 3a 8 m 2 ^d 

(n 2 - 4) 

3(n 2 - 1)K 3 _ 5a 9 m 4 ,2^ 12a 15 m 4 ,2 12 



(n 2 - 4)a 3 



K = 



-3a^ 2 - i^a 4 ; 



U = -Ka A a - a 9 m 2 (j>(j> 



(A15) 



Appendix B: Relation to the gauge invariant formalism. 



There has recently been much interest in the gauge invariant formalism 21 ' 22 which 
cast the variables of the theory (the scalar perturbations of the gravitational and scalar 
fields) into ones which are invariant under infinitesimal gauge transformations. In this ap- 
pendix we relate the different harmonics in equations (3.4) and (3.5) to the gauge invariant 



variables (in particular we shall follow the approach of [22]). 
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Mukhanov et al. define the time and space dependent scalar metric perturbations as 
ds 2 = a 2 (?7)|(l + 2<f))d7 1 2 - 2B\ i dx i dri + [(1 - 2^) 7lJ + 2E ]ij ]dx i dx j } (B.l) 
and the scalar field perturbations 

(p(x,t) = <p (t) + 6(p(x,t) (5.2) 

The above variables are related to the modes perturbations used in this paper in the 
following way 

= ^2 9nQn 



-(a n + b n )Q r 



(5-3) 

B — k n Q n 
~^(n 2 -l)V6 

~^{n 2 - 1)V6 

where we have suppressed in the sum the indices Zm corresponding to the angular momen- 



tum. 



Under a general linear gauge transformation of the form 

rj -^fj = rj + ^°(rj,x) 



(BA) 

x % —>x l = x % + 7 lJ ^|j(?7, x) 
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the scalar perturbations transform as 



a 



a' 



a 

E = E-i 



8~<j> = 8<j> - ^ ■ 



The idea of the gauge invariant formalism is to make a linear combination of the different 
scalar perturbations such that the resulting variables are independent of the gauge. A 
possible choice is 

$ = 0+ -[(B-E')aY 
a 



ct 

6<pW = dip + (p' (B — E') 
These gauge invariant quantities obey the following equations: 



Q/2 
Q/2 



which are the gauge invariant versions of the 8Gq = 8irG5T§, 8G® = &itG8Tf and 8G^ = 
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87rG5Tj equations and 

<V (si) " + mStpW _ v 2 5v? (fll) + V^aHipW - 4tf/ & + 2V^a 2 <S> = 0. (B.10) 

is the gauge invariant version of the scalar field equation. 

In the longitudinal gauge (B = k n = and E = b n = 0) used in this paper the gauge 
variables reduce to <E> = 0, \& = ip and 8(p^ 91 ' = 8(p and if we expand them in harmonics 
on the 3-sphere <E> n = g n /y/6, \l/ n = —a n /\/6 and S^ 9 ^ = fn/V®- Indeed it is easy to see 
that equations (B.9) and (B.10) are equivalent to (3.31) and (3.32) respectively, and that 
the constraint (B.8) is equivalent to (3.33). 
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